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We show that in multiband superconductors even small interband proximity effect can lead to 
a qualitative change in the interaction potential between superconducting vortices by producing 
long-range intervortex attraction. This type of vortex interaction results in unusual response to low 
magnetic fields leading to phase separation into domains of a two-component Meissner states and 
vortex droplets. 



The textbook classification of superconductors, divides 
them into two classes, according to their behavior in an 
external field. Type-I superconductors expel low mag- 
netic fields, while elevated fields produce macroscopic 
normal domains in the interior of the superconductor. 
Type-II superconductors possess stable vortex excita- 
tions which can form a vortex lattice as the energetically 
preferred state in an applied magnetic field. This picture 
of type-II superconductivity, as well as the essence of the 
more complex physics of fluctuating vortex matter, relies 
on the fact that the interaction between co-directed vor- 
tices is purely repulsive. In [1^ it was demonstrated that 
in U{1) X U{1) superconductors with two independent 
components, in a wide parameter range, there are vortex 
solutions which are on the one hand thermodynamically 
stable, and on the other hand, possess a non-monotonic 
interaction potential, repulsive at short distances but at- 
tractive at larger distances. Long range vortex attrac- 
tion in the models |l] originates from the circumstance 
that the coherence length of one of the components is the 
largest length scale of the problem, and the core of one of 
the components extends to the region where current and 
magnetic field (which are responsible for repulsive inter- 
vortex interactions) are exponentially suppressed. Indeed 
such a vortex interaction, along with their demonstrated 
thermodynamic stability, should cause the system to re- 
spond to external fields in an entirely different way from 
the vortex states of traditional type-II superconductors. 
Namely, the attraction between vortices should, at low 
fields, produce the "semi-Meissner state" [1 ) featuring 

(i) formation of voids of vortex-less states, where there 
are two well developed superconducting components and 

(ii) vortex clusters where the second component is sup- 
pressed by overlapping of vortex outer cores. This kind 
of external field-induced "phase separation" which, from 
the point of view of the second component, resembles a 
mixed state of type-I superconductors, can be interpreted 
as the system showing aspects of type-I and type-II mag- 
netic response simultaneously. The term type- 1.5 super- 
conductivity was coined for this kind of behavior. Note 
that this magnetic response originates from the existence 
of three fundamental length scales in the problem (in 



contrast to the ratio n of two fundamental length scales 
which parametrizes single-component Ginzburg-Landau 
theory), and thus it is entirely different from the inho- 
mogeneous vortex states in single-component supercon- 
ductors where inhomogeneity can be induced by defects 
in a type-II superconductor or by tiny attraction caused 
by various nonuniversal microscopic effects beyond the 
Ginzburg-Landau theory which might be pronounced in 
single-component superconductors with n is extremely 
close to l/^/2 [2]. 

Recently there has been strong and growing interest 
in multi-band materials where intercomponent interac- 
tion can be substantial. Examples are MgB2 O [4] and 
possibly new iron-based superconductors [5^. The two- 
band superconductor MgB2 [3, 4 was regarded in early 
theoretical and experimental works as a standard type-II 
superconductor. This was disputed in the recent works 
by Moshchalkov et al 0, which reported highly in- 
homogeneous vortex states formation in clean samples 
in low magnetic fields with vortex clusters (with a pre- 
ferred intervortex separation scale) and vortex-less Meiss- 
ner domains strikingly similar to the picture of the semi- 
Meissner state [1 . In connection with the experiments 
and recent suggestions that iron pnictides may also 
be multi-component superconductors, the question arises 
under what conditions type- 1.5 superconductivity is pos- 
sible (even in principle) in general multi-band systems 
with a substantial interband coupling. 

In this Letter we show that type- 1.5 behaviour can 
arise via a new mechanism in a two-band system with a 
direct coupling between the bands. The situation which 
we consider is, in a way, antipodal to that considered in 
pl : namely where only one band is truly superconducting 
while superconductivity in the other band is induced by 
the interband proximity effect. We address the properties 
of such a regime by studying the following free energy 
density (in units where ft = c = m = 1 and e is the 
Cooper pair charge). 
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Here ?/^i^2 represent the superconducting components as- 
sociated with two bands. The radical difference with pre- 
vious studies pi is that in ([T]) the effective potential for 1^)2 
has only positive terms a, /3 > 0, i.e. this band is above its 
critical temperature. It has a nonzero density of Cooper 
pairs only because of the interband tunneling represented 
by the term — r/l'^i 1 1'02 1 cos(^2 — ^1) (since the Josephson 
term favours locked phases we have 9i =62 = 0). The 
results can be generalized to including other mixed gra- 
dient and density terms in ([T]). In what follows we will 
denote the ground state values of and |V^2| by ui and 
U2. Note that in this model, in general, no explicit ex- 
pressions for ui and U2 in terms of a, /3, r] exist, but one 
can compute power series expansions for them in 7^, 
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Vortex solutions of the model take the form tpa = 
(Ta{r)e^^ ^ A = r~-^a(r)(— sin ^, cos 6>), where cro(oo) = Uq 



and a(oo) 



To understand the long-range be- 



haviour of a vortex, we choose gauge so that ipi , 1^2 are 
real, set ifji = Ui -\- Xi^ and linearize the model about 
X = {xiiX2)^ = (0,0)-^, A = 0. The result is a coupled 
Klein- Gordon system with energy density 
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where H is the Hessian of V about the ground state I'^i 
Ui, that is, Hij = d'^V/d\ilJi\d\iljj\. Clearly, 
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The eigenvalues of H are the squared masses of the nor- 
mal modes about the ground state. If 7^ = 0, then 
Xi and X2 decouple and have masses 2 and >/2a, the 
first one being in this limit the inverse coherence length 
of the first condensate, as expected. If 7^ > 0, both 
condensates have nonzero ground state values ui and 
U2 which are not known explicitly, but importantly the 
normal modes are not Xi? X2, but rather an orthogo- 
nal pair (xi cos uj + X2 sincj, — Xi sincj + X2 cos uj) where 
(coscj, sincj)-^ and (— sincj, coscj)-^ are the eigenvectors 
of H. Physically this means that the recovery of densities 
in both bands from the core singularity has a strong mu- 
tual dependence. The London penetration length is given 
by the inverse mass of A: /j^a = gy^tx^ + u^. So the linear 
theory predicts, at large r, the asymptotic formulae 

IV^il ^ ui - qi cosujKo{/iir) + q2 smujKo{/i2r) 
|V^2| ^ U2 - qi sin ujKo{jj.ir) - q2 cos ujKo{jj.2r) 
\A\ ^ r-\e-^ ^ qAKii^AT)) (5) 

where Km denotes the 7ri-th modified Bessel function of 
the second kind, and qi,q2-,qA are some unknown real 



constants depending on a,/3,r],e. Recall that Km{r) ~ 
(7r/2r)2e~^ for all m. This means that, in spite of 
the presence of two superfluid densities, we cannot talk 
about two distinct coherence lengths (in the GL sense) 
pertaining to these condensates: the leading term in both 
\^|Jl\—Ul and |'02| — '^2 decays exponentially with the same 
length scale, ^ = max{/ij"^, /i^^}. At the same time the 
system retains three fundamenatal length scales, which 
in this case are the magnetic field penetration length and 
two inverse masses /i^^, /i^^ of the modes associated with 
density variation in the coupled bands. Applying the 
methods of [8], one finds that the asymptotic interaction 
potential for two well-separated vortices is 

V ^ 2n[qlKo{^Ar) - qlM^ir) - qlKo{^2r)]. (6) 

The first term represents repulsion due to current-current 
and magnetic field interactions, while the last two terms 
represent attractive forces associated with nontrivial den- 
sity modulation, mediated in this case by the normal 
modes, described by scalar fields of mass fii and /i2- 
Hence, the linearized theory predicts that vortices should 
attract at very large separations if min{/ii,/i2} < Ma 
(and repel if min{/ii, /i2} > Ma)- It should be emphasized 
that the above analysis concerns the leading asymptotics 
of the vortex fields, not the core structure of the vortex 
directly. As we discuss below, in the presence of in- 
terband Josephson tunneling the detailed core structure 
is principally important for the form of the vortex inter- 
action potential (in contrast to usual single-component 
superconductors). This core structure cannot be de- 
rived in a simple manner from eq. ([T]). Rather, it must 
be deduced from numerical solutions of the full, nonlin- 
ear field equations which are presented in the second half 
of the paper. 

However for small rj we can demonstrate analytically 
that vortices in the model ([T]) can attract one another 
at long range but repel at short range by finding the 
masses /iA,Mi,M2- That is, we can find expansions for 
these masses and the "mixing angle" cj, valid for small 7^, 



MA = e + 0{r]^); Mi = 2 + 0{r]^); M2 



2a + 0(7?') 



00 = 



\2a-M 



(7) 



There is a range of parameters where Ui^jHi and uj can 
be computed explicitly, namely /3 = 0. Physically, this 
limit is sensible if |?/^2| remains everywhere small, which 
it does for small 7^, since U2 = 0{r]). One finds that 
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where p = 8a -h 4a^ -h 3?^^ . To understand the intervortex 
forces at short range we note that for small 7/, ip2 remains 
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FIG. 1. (Color online) Interaction energy between two vor- 
tices as a function of vortex separation in units of 10~^Ev 
where Ey is the single vortex energy for a density ratio of 
ul/ul = 0.1. 



close to zero throughout the core of ipi and in the most 
of the flux-carrying area, so one expects ?/^2 to contribute 
negligibly to the interaction energy at short range. In 
this limit one can approximate the vortex solution at this 
scale by setting ?/;2 = in F yielding a one component GL 
model with GL parameter kgl = >/2e~^, leading one to 
predict short range vortex repulsion for < e < 2 for the 
effective potential given in eq. Q. So, linear and quali- 
tative analysis suggests that the model ([T]) does possess 
type- 1.5 superconductivity at least whenever < e < 2 
and the condition for long range attraction holds, namely 
min{/i^,/i2} < e'^{u1 + 1^2)5 where /if are the eigenval- 
ues of H. In order to test this prediction, and to study 
regimes where analytic estimates cannot be made, we 
have performed numerical studies of the model ([T]) at 
various parameter values. The computation was con- 
ducted as follows: First, two phase windings were created 
around two fixed points on a numerical grid. Then, the 
free energy was minimized with respect to all degrees of 
freedom using a local relaxation method, constrained so 
that the vortex cores positions remained fixed. The pro- 
cess was then repeated for various separations, yielding 
an intervortex interaction potential. 

First let us consider the regime where the fourth or- 
der term in |?/^2| can be neglected (i.e. f3 = 0). In this 
case we conducted computations with the density ratios 
|'^2p/|'^iP being 0.1 and 0.5. The results for intervortex 
interaction energy are presented in Fig. [l][3j The com- 
puted interaction energy is given in units of 2Ey where 
Ey is the energy of an isolated single vortex. The length 
is given in units of a/2^i where ^1 is a characteristic con- 
stant (the same for all figures) defined as the coherence 
length which can be associated with this band in the limit 
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FIG. 2. (Golor online) Interaction energy between two vor- 
tices as a function of vortex separation for a density ratio of 
ul/ul = 0.5. 

of zero coupling to the second band. 

In the first case, with density ratio 0.1 we find that 
in general the density profiles of the condensates can 
be quite different, even though one of the bands has 
proximity-induced superconductivity. This can be as- 
cribed to the fact that the mixing angle uj is small (note 
that UJ ~ — lo^^/ii ) so that the subleading normal mode 

ui \2a—4:\ / 

(of mass /ii > jj.2) dominates ipi at intermediate range. 
We find that as a consequence of the disparity in the 
recovery lengths the system crosses over from type-II to 
type- 1.5 behaviour when a and rj are sufficiently low (Fig. 
[1]). The low density of means that the attractive part 
of the interaction is weak. In the curves 3 and 4, we 
find a slight long range attraction yielding a minimum 
energy at around the separation of r = 8. In the second 
case (Fig. [2|, the density ratio is increased to 0.5. The 
vortex- vortex binding energy is now much larger, and the 
minimum energy occurs at a smaller separation. Long 
range attraction occur in curves 3-5 with a maximum a 
of 0.5, in contrast to a ^ 0.1 in the previous case. 

In the third case (Fig. [3|, the electric charge has been 
increased by a factor >/2, (which is equivalent to decreas- 
ing penetration length) which decreases the magnetic re- 
pulsion between vortices. Observe the emergence of a 
new phenomenon: now the energy of an axially symmet- 
ric vortex solution with two flux quanta is smaller than 
the energy of two infinitely separated one-quanta vor- 
tices. Nonetheless, the axially-symmetric two-quantum 
vortex is not stable since the minimum energy occurs at 
nonzero vortex separation. 

Fig. [4] shows the effect of the addition of a fourth 
order term with coefficient P in the free energy of the 
proximity-induced component. It demonstrates the per- 
sistence of type- 1.5 superconductivity when fourth order 
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FIG. 3. (Color online) Interaction energy between two vor- 
tices as a function of vortex separation for a density ratio of 
u\lu\ = 0.5 for e = 1.41. 



D) 


8 


Q) 




en 


6 


X 




rte 


4 






§ 






2 
















C 











-2 


o 




o 

CO 


-4 







_g 


-6 


X 







-8 


o 




> 


-10 










X 10 



Curve 

- 1 

- 2 
3 
4 



OL (3 


V 


e 


i^2r/i^ip 




0.1 1 


1.09 


1.41 


0.5 


2.06 


0.25 0.1 


0.35 


1.41 


0.36 


1.5 


0.1 0.5 


0.55 


1.41 


0.5 


1.78 


0.1 0.1 


0.215 


1.41 


0.5 


1.61 




6 8 
Intervortex distance 



10 



12 



14 



FIG. 4. (Color online) Intervortex interaction in the presence 
of fourth-order term for ip2 in various regimes. 



terms are present for • 

To illustrate the actual behaviour of the fields lead- 
ing to this unusual intervortex interaction we plot in fig. 
[5] cross-sections of the density and magnetic field profile 
corresponding to parameter set 2 in Figj4j The figure 
clearly shows that, in spite of the identical long range 
asymptotics of density behaviour in both bands (as pre- 
dicted by the linear theory), the rate of density recovery 
in both bands at intermediate scales is actually different. 

In conclusion, we considered vortex matter in a situ- 
ation which can take place in two-band systems: only 
one band is superconducting while superfiuid density is 
induced in another band via an interband proximity ef- 
fect. This situation is in a way antipodal to the previ- 
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FIG. 5. (Color online) The behaviour of iV^il (curve 1), |'02| 
(curve 2) and magnetic field (curve 3) for a = 0.25, 77 = 
0.35, /3 = 0.1, e = 1.41. Separations are (a): 4.24 (cor- 
responding to attractive, | '02 1 -dominated interaction); (b): 
2.83 (vicinity of the minimum of the interaction potential); 
(c): 2.12 (corresponding to domination of repulsive current- 
current and magnetic interactions). 



ously studied unusual vortex interaction arising in con- 
densates with independent coherence lengths jl], e.g. as 
we showed, the asymptotics of the superfiuid densities at 
large distances from the core in both bands are governed 
by the same exponential law. However we find that, in 
contrast to the conventional single- component situation, 
the presence of even a tiny interband proximity effect can 
be crucially important. Namely, it gives rise to three 
fundamental length scales in the problem and to non- 
trivial variations of the relative superfiuid densities in 
two bands in a vortex producing type- 1.5 behaviour in a 
wide range of parameters. It should manifest itself in the 
magnetic response which involves a phase separation into 
vortex and two-component Meissner domains. The effect 
may be more common near the temperature where the 
weak band crosses over from active to proximity-induced 
superconductivity because a should be small near this 
temperature. 

We thank Alex Gurevich, V. Moshchalkov and M. 
Wallin for discussions. The work is supported by the 
Swedish Research Council, NSF CAREER Grant No. 
DMR-0955902, and the U.K. Engineering and Physi- 
cal Sciences Research Council. E. B. was supported by 
the Knut and Alice Wallenberg Foundation through the 
Royal Swedish Academy of Sciences. 



[1] E. Babaev & J.M. Speight Phys.Rev. B 72 180502 (2005) 
[2] E. H. Brandt, Rep. Prog. Phys. 58, 1465 (1995) 
[3] A. Liu, I.I. Mazin, J. Kortus, Phys. Rev. Lett. 87 087005 
(2001); 

[4] A. Gurevich, Phys. Rev. B 67 184515 (2003) ; Physica C 
056160 (2007) 

[5] see e.g. K. Ishida, Y. Nakai,H., Hosono, J. Phys. Soc. Jpn. 

78 062001 (2009) 
[6] V.V. Moshchalkov, al Phys. Rev. Lett. 102, 117001 

(2009) 



5 



[7] T. Nishio et al Phys. Rev. B 81, 020506(R) (2010) 
[8] J.M. Speight, Phys. Rev. D 55 3830 (1997) 



